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ABSTRACT 

We introduce a generalized scaling law, A/tot = 10 K A a B b , to look for the minimum scatter 
in reconstructing the total mass of hydrodynamically simulated X-ray galaxy clusters, given 
gas mass Af gas , luminosity L and temperature T. We find a locus in the plane of the loga- 
rithmic slopes a and b of the scaling relations where the scatter in mass is minimized. This 
locus corresponds to bm = — 3/2<zm + 3/2 and b^ = — 2a^ + 3/2 for A = M gas and L, 
respectively, and B = T. Along these axes, all the known scaling relations can be identified 
(at different levels of scatter), plus a new one defined as M to t oc (LT) 1 / 2 . Simple formula to 
evaluate the expected evolution with redshift in the self-similar scenario are provided. In this 
scenario, no evolution of the scaling relations is predicted for the cases (pM = 0, aj; = 1) 
and (pL — 7/2, o,l = —1), respectively. Once the single quantities are normalized to the av- 
erage values of the sample under considerations, the normalizations K corresponding to the 
region with minimum scatter are very close to zero. The combination of these relations allows 
to reduce the number of free parameters of the fitting function that relates X-ray observables 
to the total mass and includes the self-similar redshift evolution. 

Key words: cosmology: miscellaneous - galaxies: clusters: general - X-ray: galaxies: clus- 
ters. 



1 INTRODUCTION 

Galaxy clusters are believed to form under the action of gravity in 
the hierarchical scenario of cosmic structure formation (e.g. Voit 
2005). They assemble cosmic baryons from the field and heat them 
up through adiabatic compression and shocks that take place dur- 
ing the dark matter halo collapse and accretion. Simple self-similar 
relations between the physical properties in clusters are then pre- 
dicted (e.g. Kaiser 1986, 1991, Evrard & Henry 1991) since gravity 
does not have any preferred scale and hydrostatic equilibrium be- 
tween intra-cluster medium (ICM) emitting in the X-rays (mostly 
by thermal bremsstrahlung) and the cluster potential is a reason- 
able assumption. These scaling relations are particularly relevant 
to connect observed quantities, such as X-ray luminosity, temper- 
ature and mass, to total cluster mass, which is used to constrain 
cosmological parameters (e.g. Allen, Mantz & Evrard 201 1). 

Work in recent years has focused in defining X-ray mass prox- 



ies, i.e. observables which are at the same time relatively easy to 
measure and tightly related to total cluster mass by scaling rela- 
tions having low intrinsic scatter as well as a robustly predicted 
slope and redshift evolution (e.g. Kravtsov et al. 2006, Maughan 
2007, Pratt et al. 2009, Stanek et al. 2010, Short et al. 2010, Fab- 
jan et al. 2011). An important role in defining such proxies and 
assessing their robustness is played currently by cosmological hy- 
drodynamical simulations, thanks to their ever improving numeri- 
cal resolution and sophistication in the description of the physical 
processes determining the ICM evolution (e.g. Borgani & Kravtsov 
2009). 

In this letter, we present and discuss the behaviour of the scal- 
ing relations generalized to include the dependence upon two inde- 
pendent observables, one accounting for the gas density distribution 
(namely gas mass A/ gas and X-ray luminosity L), the other tracing 
the ICM temperature, T. This paper is organized as follows. In Sec- 
tion 2 we introduce the scaling relations investigated. In Section 3, 
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we discuss the redshift evolution and the normalization of these 
relations, and how they depend on the selection adopted to define 
the sample analyzed. In Section 4, we summarize and discuss our 
results in view of their application to observational data. 



2 THE GENERALIZED SCALING LAWS 

Under the assumptions that the smooth and spherically symmetric 
intra-cluster medium (ICM) emits by thermal bremsstrahlung and 
is in hydrostatic equilibrium with the underlying gravitational po- 
tential, the self-similar (SS) scenario relates bolometric luminosity, 
L, gas temperature, T, gas mass, M gas , to the total mass, M to t in 
a simple and straightforward way. For instance, the equation of hy- 
drostatic equilibrium, d(p gas T) /dr « p gas GMtot/r 2 , allows to 
write Mtot oc TR, as long as the slope of temperature and gas 
density profiles are independent of cluster mass. By combining it 
with the definition of the total mass within a given overdensity 
A z with respect to the critical density at the cluster's redshift z, 
Mtot oc E z A z R' i , one obtains that E z Al /2 M to t oc T 3/2 , where 
E z = H z / H = [fi m (l + zf + 1 - fi m ] 1/2 for a flat cosmol- 
ogy with matter density parameter fj m , cosmological constant and 
Hubble constant at the present time Ho. Similarly, the definition 
of the bremsstrahlung emissivity e oc A(T)n 2 as oc T 1 ' 2 ng as (the 
latter being valid for systems sufficiently hot, e.g. > 2 keV) relates 
the bolometric luminosity, L, and the gas temperature, T: L « eR 3 
pa T 1/2 / 2 as M 2 ot i?~ 3 * /g as T 2 , where we have made use of the 
above relation between total mass and temperature. 

By combining these basic equations, we obtain that the scaling 
relations among the X-ray properties and the total mass are (see 
also Ettori et al. 2004): 

• E z M t ot oc T 3/2 

• Mtot oc M gas 

. E z Mtot oc (E- 1 i) 3/4 . 

Kravtsov et al. (2006) introduced the Yx mass proxy, which 
is given by the product of temperature and gas mass. Owing to its 
definition, it is related to the total thermal energy of the ICM. They 
demonstrated that, among the known mass indicators, Yx is a very 
robust mass proxy. Its scaling relation with M500 being charac- 
terized by an intrinsic scatter of only 5-7 per cent at fixed Yx, 
regardless of the dynamical state of the cluster and redshift, with a 
redshift evolution very close to the prediction of self-similar model. 
Arnaud et al. (2007) used XMM-Newton data of a sample of 10 re- 
laxed nearby clusters spanning a Yx range of 10 13 — 1O 15 M keV, 
and confirmed that the M500 — Yx relation has a slope close to the 
self-similar value of 3/5, independent of the mass range consid- 
ered. They showed that the normalisation of this relation is about 
20 per cent below the prediction of numerical simulations which 
include cooling and supernova (SN) feedback, and explained this 
offset with two different effects: an underestimate of true mass 
due to a violation of the assumption of hydrostatic equilibrium, 
and an underestimate of hot gas mass fraction in the simulations 
(see also Zhang et al. 2008). They confirmed that Yx might in- 
deed be a better mass proxy than T and M gas by comparing the 
functional form and scatter of the relations between different ob- 
servables and mass. Extensive use of the Yx — M to t relation has 
been made in recent analyses aimed at constraining cosmological 
parameters through the evolution of the cluster mass function (e.g. 
Vikhlinin et al. 2009) and the properties of the scaling relations 
(Mantz et al. 2010). Pratt et al. (2009) presented the X-ray luminos- 
ity scaling relations of 3 1 nearby clusters from the Representative 



XMM-Newton Cluster Structure Survey (REXCESS), all having 
temperature in the range 2-9 keV and selected in X-ray luminos- 
ity so as to properly sample the cluster luminosity function. Their 
analysis showed that scaling relations between bolometric X-ray 
luminosity and temperature, Yx and total mass, are all well repre- 
sented by power-law shapes with slopes significantly steeper than 
self-similar predictions. They concluded that structural variations 
have little effect on the steepening, whereas it is largely affected 
by a systematic variation of the gas content with mass. Maughan 
(2007) analysed Chandra ACIS-I data for 115 galaxy clusters at 
0.1 < z < 1.3 observed to investigate the relation between lumi- 
nosity and Yx • They found that the scatter is dominated by cluster 
cores, and a tight Lx — Yx relation (11 per cent intrinsic scatter 
in Lx) is recovered if sufficiently large core regions (0.15-Rsoo) 
are excluded. The tight correlation between Yx and mass and the 
self-similar evolution of that scaling relation out to z = 0.6 is con- 
firmed. Fabjan et al. (201 1) analysed an extended set of cosmolog- 
ical simulations of galaxy clusters, and confirmed that the M — Yx 
scaling law is the least sensitive to variations of the physics in the 
ICM and very close, in terms of slope and evolution, to predictions 
of the self-similar model. They also pointed out that M — M gas is 
the relation with the smallest scatter in mass, whereas M — T is the 
one with the largest among the considered scaling relations. 

In the present work, we generalise the definition of the Yx 
mass proxy, by considering the scaling relation between total mass, 
Mtot, and a more general proxy defined in such a way that M to t oc 
A a B b , where A is either M gas or L and B = T. The use of this 
relation generalizes the relation M to t — Y, while maintaining the 
attitude to recover total mass by combining information on depth 
of the halo gravitational potential (through the gas temperature T) 
and distribution of gas density (traced by M gas and X-ray luminos- 
ity), the latter being more affected by the physical processes deter- 
mining the ICM properties. In doing that, we aim to minimize the 
scatter in the relations between total mass and observables by (i) 
relaxing the assumption of the self-similarity, (ii) adopting a gen- 
eral and flexible function with a minimal set of free parameters, 
(iii) offering a method that can be readjusted in dependence of the 
specific sample selection adopted. 

In the recent past, similar work has been done by different au- 
thors with the aim of generalising the use of simple power-law scal- 
ing relations between cluster observables and total mass. Stanek et 
al. (2010) discussed the second moment of the halo scaling rela- 
tions by investigating the signal covariance at fixed mass in numer- 
ical simulations. Okabe et al. (2010) used a small sample of 12 ob- 
jects observed with Subaru and XMM-Newton to study the covari- 
ance between the intrinsic scatter in M to tMg as and M to tT rela- 
tions and to propose a method to identify a robust mass proxy based 
on principal component analysis. Rozo et al. (2010) presented an 
extensive discussion on the relaxation of some assumptions on the 
parametrization of the relation between optical richness and total 
mass, by introducing the possibility of deviation from a power-law 
shape, as well as richness- and mass-dependence of intrinsic scat- 
ter. 

To study the behaviour of the Mtot oc A a B b relation in min- 
imizing the scatter, we used a sample of 24 Lagrangian regions, se- 
lected around the most massive clusters with a radius equal to five 
times the virial radius, and extracted from a parent low-resolution 
N-body cosmological simulation with a box of size 1 /t _1 Gpc co- 
moving, as described in Bonafede et al. (201 1). A flat ACDM cos- 
mological model with Q, m = 0.24, £lt> ar = 0.04, n s — 0.96, 
as = 0.8 and present day Hubble constant of 72 km s _1 Mpc -1 , 
consistent with WMAP-7 cosmological parameters (Komatsu et 
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Figure 1. Contour plots that enclose 1.2, 1.5 and 2 times the minimum scatter, as function of the slopes a and b of the generalized scaling relations, as indicated 
in each panel. Also overplotted are the lines from eq.\T\(green dashed) and eq.[5](fc/ue dotted). (Top panels) The case of {A = Af gas , B = T} using (from 
left to right) observational-like measurements of the xmas csf sample; direct measurements of the csf sample; direct measurements of the nr sample. (Bottom 
panels) The same as above, but for the case of {A = L, B = T}. The insets show the values of the normalization if as a function of the slope a in the region 
enclosed within 2 times the minimum scatter. 



al. 2011), was assumed. A set of 24 Lagrangian regions, centred 
around as many massive clusters, were re-simulated by increas- 
ing mass resolution and adding high-frequency modes to the power 
spectrum (Tormen et al. 1997). Within the high resolution region, 
dark matter particles have a mass uibm = 8.47 x 10 s h' 1 M Q . 
The size of each Lagrangian region was chosen in such a way that 
by z = there are no low-resolution particles within at least 5 
virial radii from the central cluster. As a result, the large extent of 
each of these high-resolution regions allows one to identify more 
than one single cluster-sized halo within it, which is not contami- 
nated by low-resolution particles within its virial region (Bonafede 
et al. 2011; Fabjan et al. 2011). 

Clusters identified from this set of initial conditions were sim- 
ulated with the TreePM-SPH GADGET-3 code, an improved ver- 
sion of the original GADGET-2 code (Springel 2005). As described 
by Fabjan et al. (2011), simulations have been carried out for two 
different prescriptions for the physics determining the evolution 
of cosmic baryons: (i-sample nr) non-radiative physics and (ii- 
sample csf) including metallicity-dependent radiative cooling, a 
model for star formation and galactic winds triggered by SN explo- 
sions (as described by Springel & Hernquist 2003) with velocity 
v w — 500 km s _1 , and a detailed model of chemical evolution as 
described by Tornatore et al. (2007). By selecting only objects with 
mass weighted temperature T > 2 keV, we end up with 41 ob- 



jects in each sample. A subset of the c.s/sample has been processed 
through the X-MAS tool (e.g. Rasia et al. 2008) to generate Chan- 
dra mock observations, and then analyzed with an observational- 
like approach to measure temperatures and gas masses (xmas sam- 
ple; Rasia et al. 2011). The latter sample includes all the clus- 
ters with spectroscopic-like temperature larger than 2 keV, and ob- 
served along 3 orthogonal projection directions, so that we end up 
with 159 mock observations of simulated clusters. Total and gas 
masses within -R500 are computed as as described in Fabjan et al. 
(201 1) and Rasia et al. (201 1). Gas temperatures and luminosities, 
both bolometric and in the 0.1-2.4 keV band, are computed after 
excising cluster core regions, which are defined as the regions en- 
closed within 0.15i?5oo. The effect of core excision is also consid- 
ered in the discussion of the results and is shown not to affect the 
conclusions of our analysis. 

We fit a linear relation to the log-log scaling between to- 
tal mass and proxies, normalized to the average values computed 
within each sample of simulated clusters: 

k>gio A ^tot — K + a log 10 A + b log 10 B, . (1) 

Here we defined {iWtot = M tot /M tot ,A = A/A,B = B/B}, 
with barred quantities indicating the average values of the corre- 
sponding quantities. 

Within each set of simulated clusters, containing N ob- 
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Figure 2. Evolution of the scaling relation parametrized through the quantity E z in the SS scenario and including the relations in equation [2] Solid (dashed) 
lines show the behaviour for the relations with A = M gas (A = L). (Left) Values of c, exponent of E z , from equation|4]as function of the logarithmic slope 
b (c = in the case of no-evolution). (Right) Values of E z as function of the redshift for different scaling relations. As representative cases, two sets of lines 
are plotted: thin lines assume b = — 1, thick lines adopt 6 = 3. 



jects, we compute for each pair of values of the slopes {a,i,bj} 
the corresponding scatter, which is defined as a 2 (ai,bj) = 
X^=i,iv( lo gio^tot,k - K - a 4 log 10 i fe - bjlog 10 B k ) 2 /N, 

where K = Z) fc =i,iv( lo gio A ^tot,k - a,i\og 10 A h - 
bj log 10 Bk)/N. We then search find the locus in the {a, b} 
plane where scatter is minimized in a similar. In all cases, this 
locus is well represented by the lines 

{A = A/ gas , B = T} => 6 M = -3/2a M + 3/2 

{A = L, B = T} => b L = -2a L + 3/2, (2) 

(see Fig.QJ or, in a more concise form, 6 = — (1 + l/2d) a + 3/2, 
where d corresponds to the power to which the gas density appears 
in the formula of the gas mass (d = 1) and luminosity (d = 2). 
This correlation between logarithmic slopes allows us to reduce by 
one the number of free parameter in the linear fit of the generalized 
scaling law between observables and total mass. 

It is worth noticing that these relations reduce to the standard 
self-similar predictions in the appropriate cases: M t ot oc T 3/2 , 
Mtot oc Mg as , Mtot oc V 3 / 5 are recovered for ajvf = 0, 1 and 3/5, 
respectively; Mtot oc L 3/4 and M tot oc (LT) 1/2 , which is the 
corresponding relation of Mtot oc y 3//5 once gas mass is replaced 
by luminosity, are recovered for a_L =3/4 and 1/2, respectively. 

However, to represent the tilted shape of the contours encir- 
cling the region with the minimum scatter in the simulated dataset 
here investigated, we should prefer the following relations among 
the logarithmic slopes, 

b M ~ -1.90M + 1-8, b L w ~2Aa L + 1.8, 

6» -(1.4 + 0.5 d) a + 1.8 (3) 

that are shown as dotted lines in Fig.Q] 

In the following discussion, we refer to the SS case described 
from the equations|2]as the reference one. 



3 EVOLUTION, NORMALIZATION AND ROBUSTNESS 
OF THE GENERALIZED SCALING LAWS 

In this section, we discuss some properties on the redshift evolution 
and normalization of the generalized scaling laws, and present the 
results of the tests by which we have verified the robustness of our 
predictions. 

3.1 Evolution of the generalized scaling laws 

With simple mathematical substitutions, we can predict the redshift 
evolution expected for the SS case, Aftot oc E c z : 

{A = M gas , B = T} => c M = -2/36 M = a M - 1 

{A = L, B = T} => c L = b L /2 - 7/4 = -a L - 1. (4) 

We are now in the position to look for the scaling relation 
which has the weakest redshift dependence or, on the contrary, the 
relation which makes this dependence stronger. We note that there 
is no dependence on redshift only in two cases among the scaling 
relations here investigated (see Fig.|2](: (i) om = 1 (and bM = 0), 
i.e. for the scaling law Mtot oc M gas ; (ii) a_L = — 1 (and 6_l = 
7/2), i.e. for the relation M to t oc Z/ -1 T 7//2 . The prediction for the 
lack of evolution of these scaling relations can be tested against 
observational data. 

3.2 Normalization of the generalized scaling laws 

As shown in Fig.[T] the normalization K corresponding to the value 
of minimum scatter is close to zero. This is expected once the quan- 
tities are normalized to the averaged values Mtot, A, B. However, 
only A and B are known for an observed sample. Thus, by adopt- 
ing one of the relations in equation [2] one can directly measure 
Mot =■ Mtot /Mtot and recover the total mass Mtot only once 
Mtot is independently evaluated either through mock samples se- 
lected from catalogs of hydrodynamically simulated objects to con- 
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Figure 3. Best-fit values of the slopes clm and 6m as a function of the sample examined, (1)-(10) from xmas (Rasia et al. 2011), (11)-(14) from direct 
measurements in the hydrodynamical simulations (Fabjan et al. 201 1): (1) at A = 2500, with the spectroscopic-like estimate T = T s \ > 2 keV; (2) at 
A = 2500, with the X-ray spectroscopically determined T = T x > 2 keV; (3) with T = T x > 2 keV and the core included; (4) with T = T s \ > 2 keV; 
(5) with T = T x > 2 keV; (6) with T = T sl > 4 keV; (7) with T = T x > 4 keV; (8) with A/tot > 10 14 M Q ; (9) with T = T x and 1,000 realizations 
of randomly selected objects with 30 clusters with 2 < T < 4 keV and 40 clusters with 4 < T < 10 keV; (10) the same as in (9) but including a relative 
statistical error of 20 per cent on the total mass; (1 1) for the sample nr and excluding the core (0 — 0.15-Rsoo): (12) including the core, (13) for the sample 
c.s/and excluding the core, (14) including the core. For the cases (3)-(14), all quantities are estimated at A = 500. The bottom panels show the corresponding 
scatter a. Green diamonds and blue squares in the central (lower) panels are the predicted slope 6m (scatter) from the SS relations in equation [5] and the 
corrected relations in eq.fj] respectively. 



tain the same number of objects, and with similar properties, of the 
observed ones, or through a self-calibration tuned by a sub-sample 
of clusters for which robust mass estimates are available. Under this 
respect, the suggested approach is the standard one, with the same 
limitations affecting any other application of the scaling laws: mass 
calibration and selection effects. The innovation, we are propos- 
ing, is to add an extra parameter, imposing a new constraint on the 
slopes of the scaling laws, to allow a further minimization of the 
scatter. 



3.3 Robustness of the generalized scaling laws 

To assess the robustness of the analysis of the simulated dataset, we 
have repeated our calculations by extracting the simulated objects 
according to different criteria, e.g., including or excluding the clus- 
ter core emission, adopting different overdensity, using different 
definition for the gas temperature, selecting only very hot or mas- 
sive systems. All these samples reproduce consistently the plots 
shown in Fig. Q] by varying only the location of the best-fit val- 
ues, but confirming the dependence among the logarithmic slopes 
over the region of the parameter space that minimize the measured 
scatter (see Fig.[3]l- 

When observational data are considered, several other selec- 
tion effects can still affect both the definition of a sample and the 
measurements of the normalization and slope of the adopted scaling 
law. A proper treatment of the second-order moments and of the 
covariance related to the scaling relation has then to be addressed 
(see, e.g., Stanek et al. 2010, Rozo et al. 2009 and 2010, Mantz et 
al. 2010). 



4 SUMMARY AND DISCUSSION 

We have presented new generalized scaling relations with the 
prospective to reduce further the scatter between mass proxies and 
total cluster mass. We find a locus of minimum scatter that relates 
the logarithmic slopes of the two independent variables considered 
in the present work, namely temperature T, which traces the depth 
of the cluster potential, and another one accounting for the gas den- 
sity distribution, such as gas mass M gas or X-ray luminosity L. 
Within this approach, all the known scaling relations appear as par- 
ticular realizations of generalized scaling relations. For instance, 
we introduced the scaling relation Mtot oc (LT) 1 / 2 , which is anal- 
ogous to the Mtot — Y relation, once luminosity is used instead of 
gas mass. 

Also the evolution expected in the framework of the self- 
similar model are predicted for the generalized scaling relations. 
They can be used either to maximize the evolutionary effect to test 
predictions of the self-similar models itself or, on the contrary, to 
minimize them in case of cosmological applications. 

A linear function in the logarithmic space can be then fitted to 
the data normalized to the average values measured in the sample: 

log 10 Mtot = K + olog 10 A + 61og 10 B + clog 10 E z , (5) 

with K = 0, b M = -3/2om + 3/2, c M = -2/3b M = au - 1 
for {A = Af gas , B = T} and K = 0, b L = -2a L + 3/2, c L = 
6l/2 — 7/4 = — <zl — 1 for {A = L, B = T}. In a more concise 
form, the above relation can be recast as 6 = — (1 + 1/2 d) a + 3/2, 
where d corresponds to the power with which gas density appears 
to define either gas mass (d — 1) or luminosity (d — 2). This 
fitting function has 4 free parameters that are reduced to one (plus 
the average value of the total mass of the objects in the sample) 
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thanks to the existing tight correlation found between a and 6, at 
least within the region of the {a, b} parameter space where intrinsic 
scatter is minimised. 

The method and the results presented in this work offer a ro- 
bust framework to relate, with the request of a minimum scatter, the 
X-ray observables to the total gravitational mass of galaxy clusters 
for studies of their thermodynmical properties and for cosmological 
application. 
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